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Theta vacuum effects on the pseudoscalar condensates and the rf meson 
in 2-dimensional lattice QED * 

Hidenori Fukaya a , Tetsuya Onogi ^ 

a Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

We study the chiral condensates and the rf meson correlators of the massive Schwinger model in 9 ^ vacuum. 
Our data suggest that the pseudoscalar operator does condense in a fixed topological sector and gives long range 
correlations of the rf meson. We find that this is well understood from the clustering decomposition and statistical 
picture. Our result also indicates that even in 8 = case, the long range correlation of rf meson receives non- 
zero contributions from all the topological sectors and that their cancellation is non-trivial and requires accurate 
measurement of the reweighting factors as well as the expectation values. It is then clear that the fluctuation of 
the "disconnected" diagram originates from the pseudoscalar condensates. 



1. Introduction 

In QCD or the massive Schwinger model in 
6 = vacuum [1-4], it is well known that the 
scalar operator condense while the pseudoscalar 
does not; 



(M^O, (^75 V) =0 {0 = 0), 



(1) 



where the second equation follows from parity 
symmetry. However, if we have a non-zero 6 
term, which violates parity symmetry, both of 
them have non-zero expectation values; 



(W) ? 0, (^75 V-) ^0 (0^0), 



(2) 



which indicates that rj meson should have a long- 
range correlation as 

lim (vHxHO)) cx (^75^(0:)) (^75^(0)) ^ 0. (3) 

| as | — >large 

We would like to present our numerical results 
of the 2-flavor massive Schwinger model with a 
0^0 term. We investigate tp^ip condensates 
and the rf meson correlators in each topological 
sector. It is also found that they are non-trivially 
related each other to reproduce the 9 dependence. 
We find that their behavior is well understood 
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by the intuitive picture based on the clustering 
decomposition and the statistical mechanics. 

In particular, our study shows that the accu- 
rate contributions from higher topological sectors 
are essential in order to assure parity symmetry, 
which never allows the long-range correlation of 
the rf meson correlators. Moreover, it is also 
shown that the origin of the fluctuations of dis- 
connected diagram is from pseudoscalar conden- 
sates in each topological sector. 

2. Strategy for simulations 

Our strategy to calculate the 9 vacuum effects 
is to separate the integral of the gauge fields into 
topological sectors; 



(O) 



f3,m 



(4) 



where f3 = 1/g 2 denotes a coupling constant and 
(0)Q and denote the expectation value and 
the partition function in a fixed topological sector 
respectively. 

The expectation values with a fixed topologi- 
cal charge; (O)® are evaluated by generating link 
variables with the following gauge action [5]; 



Sg = 



^l-(l-ReP^ (*))/£ 

if admissible 



3C 



otherwise, 



2 



where P M „ denotes a plaquette and e is a fixed 
constant. This action impose the Luscher's 
bound [5] on the gauge fields, which realize an ex- 
act topological charge on the lattice, that is never 
changed in each step of the the hybrid Monte 
Carlo updation. 

2jQ normalized by that of zero topological sec- 
tor can be evaluated by decomposing it into three 
parts; 



[P ' m) ~ Z°(J3,m) * ' 

classical solution 

J cL4g det(£> + m) 2 J~ dp' AS° ,m) 
X JdA° cl dct{D + m) 2 ^ 



(6) 



fluctuation 



moduli integral 



where SQ min denotes the classical minimum of 
the gauge action with topological charge Q, 
J dA^ denotes the moduli integral, and AS® = 
( S g - S2 min )% - (S G )% , all of which arc numer- 
ically calculable [6,?]. 

We choose the domain- wall fermion action with 
Pauli-Villars regulators for sea quarks. The link 
variables are updated by the hybrid Monte Carlo 
algorithm. The parameters are chosen as = 1.0, 
m = 0.1,0.15,0.2,0.25,0.3. We take L x = L y = 
16 and L s = 6 lattice where L s denotes the size 
of the extra dimension of domain-wall fermions. 
50 molecular dynamics steps with a step size 
At = 0.035 are performed in one trajectory. Con- 
figurations are updated every 10 trajectories. For 
each topological sector, around 500 trajectories 
are taken for the thermalization staring from the 
initial configuration which is the classical instan- 
ton solution with topological charge Q. We gener- 
ate 300 configurations in —\Q m ax\ <Q< \Qmax\ 
sectors for the measurements and from 1000 to 
10000 for the rewcighting factors at various 0, 
where |<3 ma x| = 

4 at m = 0.1,0.15,0.2 and 

\Qmax\ — 5 at TTL — 0.25,0.3. 

3. Results 

The topological charge dependence of pseu- 
doscalar condensates is derived by the anomaly 
equation; 

-(^} Q =Q/(mV), 



where V denotes the volume of the torus. As seen 
in Fig. 1, our data show a good agreement with 
this equation. 

Then the 77' meson correlators should have 
long-range correlations. From the clustering de- 
composition, this can be expressed as 

- E p Q>Q' <E ^n^f) Q B <E ^n^f)T Q ' . (8) 

Q' f 1 

where ()^ B means the expectation value with the 
topological charge Q' in the region A, B which 
denote the half of the large box, where the pseu- 
doscalar operators reside, respectively and Pq,q> 
denotes the probability of the distribution where 
Q' instantons appear in the box B and Q — Q' 
appear in the box A. 

In Q = case, one obtains 



\v ( x )v (0)) -> 

- E P °>Q' E i>n^f)% (E ^n^sY 

Q' f f 



Q' 



= + E p w ^E ^n^ffs j > 0. ( 9 ) 

where we assume (0)a — (O)b and use the anti- 
symmetry; (^75 ■0) _< 3 = —(ipj5ip)Q as seen in 
Fig. 1. On the other hand, at large Q, assum- 
ing the distribution Pq,q> to be Gaussian around 
Q' <~ Q/2, the correlation can be evaluated as 
follows, 



Q |a;|^largo 



-V^A J dQ'e- a W-W 2 



x (E^/^ 5 ^/)b (E^/^ 5 ^/)a Q ' 
/ / 



i 2 V 2 



iQ'/mV ~4,(Q-Q')/mV 

4Q 2 



4 / n2 1 



a I m 2 V 2 



(10) 



(7) 



where 1/y/a « Q/2 is a numerical constant. 
As seen in Fig 2, it is surprising that these very 
simple arguments describe the data quite well. 
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9 dependence of the rj' correlators are evaluated 
by substituting the data into Eq.(4). Fig. 3 shows 
the result. It is obvious that there are long-range 
correlations at 9 ^ while 9 = case is consistent 
with zero, which suggests our reweighting method 
works well at small 9. 

4. Summary 

We study Q and 9 dependence of the pseu- 
doscalar condensates and the rj' meson correla- 
tors. We find that pseudoscalar does condense in 
each topological sector; 



(ii) 



and there exists a long-range correlation of rj' me- 
son; 

lim { ?7 ' t (a;)r/'(0)) Q=o > 0, 

|x|— s-largc 

lim (flt( I W(0)p 2 ~-^<0 ) (12) 

|x|^large TJl^V z 

which are well understood by the clustering prop- 
erties. 

It is also found that each contribution from 
different topological sectors plays very important 
role to produce non-trivial 9 dependence of these 
observables. In particular, the cancellation the 
long-range correlation of rj' meson requires accu- 
rate measurements of higher topological sectors. 
It is also obvious that the fluctuation of the dis- 
connected diagrams originates from these pseu- 
doscalar condensates. 

It would be interesting to extend our studies to 
4-dimensional QCD. 
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Figure 1. — (^75^) Q in each topological sector 
at (3 = 1.0. The lines show the anomaly equa- 
tion — ("ip j5ip}® = Q/(mV), which agrees with 
our data. 
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Figure 2. The long-range correlations in each 
topological sector evaluated by the fit with con- 
stants. The lines are (r/t(x)?/(0)) Q l^ argc 
-4Q 2 /(mV) 2 . 
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Figure 3. The propagations of the 7/ meson at 
6 = 0,0. 37T and m = 0.1 are shown. The lines 
show the result of the fit. 



